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We have considered the S — 1/2 antiferromagnetic Heisenberg model in two dimensions, with 
an additional Ising next-nearest neighbour interaction. Antiferromagnetic next-nearest neighbour 
interactions will lead to frustration, and the system responds with flipping the spins down in the xy 
plane. For large next nearest neighbour coupling the system will order in a striped phase along the z 
axis, this phase is reached through a first order transition. We have considered two generalizations 
of this model, one with random next-nearest neighbour interactions, and one with an enlarged 
unit cell, where only half of the atoms have next-nearest neighbour interactions. In both cases 
the transition is softened to a second order transition separating two ordered states. In the latter 
case we have estimated the quantum critical exponent j3 ~ 0.25. These two cases then represent 
candidate examples of deconfined quantum criticality. 

PACS numbers: 75.10.Jm,75.10.Nr,75.40.Mg,75.40.Cx 



I. INTRODUCTION 

The groundstate of the antiferromagnetic Heisenberg 
model is macroscopically degenerate. This makes it par- 
ticularly sensitive to additional interactions, which might 
induce transitions to different states In a seminal pa- 
per the concept of fractionalizcd order, was set forth in 
Ref. Q • The generic starting point of this analysis is the 
two dimensional antiferromagnetic Heisenberg model 

H = jJ2Si-Sj + --- , (1) 

(id) 

where the ellipsis represent additional short range inter- 
actions, governed by a coupling g. For g = the ground- 
state is the antiferromagnetic Neel state, and by tuning 
g the system can supposedly be driven through a contin- 
uous quantum phase transition to a state with a different 
type of order. According to the Landau-Ginzburg- Wilson 
(LGW) paradigm for phase transitions, an order-to-order 
transition must either be first order, or proceed via an in- 
termediate disordered state. Since the scenario envisaged 
in Ref. 0, namely a continuous order-to-order transi- 
tion, breaks with this paradigm, the term "deconfined 
criticality" was coined to describe these transitions ■ 

The only microscopic model considered in some detail 
in the context of deconfined criticality is a dimer model 
with two spins in the unit cell[j|. For this particular 
model the transition between a Neel state and a spin- 
gapped paramagnet can be shown analytically, and are 
also confirmed with QMC calculations 3 . Apart from 
this dimer model it has been difficult construct micro- 
scopic models which give rise to deconfined criticality. 
Sandvik et.al. have investigated a model with ferromag- 
netic XY interactions and a four-spin ring exhange0,|g, 
this model has a quantum critical point separating a su- 
pcrfluid and valence bond solid, which might be a micro- 
scopic manifestation of deconfined criticality. Another 
possible way to build a microscopic model which might 



give rise to deconfined criticality is to include frustration. 
A natural way to frustrate the Heisenberg model is with 
a next nearest neighbour (nnn) Heisenberg interaction; 
this is usually called the Ji — J2 model. The J2 coupling 
will favor antiparallel spins along next nearest neighbour 
bonds, this is in conflict with the nearest neighbour ex- 
change. The result is frustration, and a reduction in the 
antiferromagnetic ordering. 

Unfortunately, in this model the geometric frustra- 
tion gives rise to a sign-problem, and the model is re- 
ally not amenable to a Monte Carlo based approach. 
Studies of this model have been based on a reweight- 
ing technique^, exact diagonalizationQ and variational 
methods 9] . The results indicate that Neel order persists 
up to k — J2/J1 < 0.40, and that a striped order develops 
for k > 0.60. Recent results indicate that the transition 
at k ~ 0.40 is a weak first-order transition 10] . 

To avoid the sign problem of the J\ — J2 model, one 
can study a simplified model where the nnn exchange is 
only along the z-components of the spin, i.e. the model 

H = J^Si-gj + K J2 s i s !)- W 

(id) ((id)) 

This simplified model captures the effect of frustration, 
but in contrast to the J\ — J2 model the isotropy in spin 
space is explicitly broken by the next-nearest neighbour 
interaction. Hence, in particular for small k we expect 
this to be a much stronger perturbation of the antiferro- 
magnet Heisenberg model than the J2 coupling. 

Apart from the Heisenberg point at n — we expect 
three different phases as k is varied: For k < the sys- 
tem is not frustrated, and the additional next-nearest 
neighbour will only serve to increase the antiferromag- 
netic ordering. Observe however that the next-nearest 
neighbour interaction has singled out the z direction in 
spin space, i.e. the model should be in the universal- 
ity class of the Ising model and have an ordered state 
at finite T. For k > the system will be frustrated, for 



2 



moderate k we expect that the system will avoid the frus- 
tration by flipping the spins down in the xy plane, i.e. 
we will effectively get an antiferromagnetic 0(2) model. 
For larger values of k the next-nearest neighbour inter- 
action will dominate, in which case the spins will again 
point along the z axis, and order in a state with stripe 
order. The transition from the effective antiferromag- 
netic 0(2) model to the striped state is first order. In 
Ref. 1 1 11 ] an antiferromagnetic Heisenberg model with an 
additional anisotropic next-nearest neighbour exchange 
was studied. This work reported Monte Carlo results in 
the limit of zero transverse nnn interactions; i.e. Eq. [21 
In addition, a model similar to Eq.[2]in the terms of hard- 
core bosons was considered in Refs. E3 We have 
deteremined the phase diagram of Eq. essentially re- 
producing the results of Refs. |llj, |l2j, |l3|. In addition, we 
have considered two generalisations aimed at softening 
the first order transition to the striped state. 

The generalisations of Eq. [21 we have considered are 
first a disordered model, where the next-nearest neigh- 
bour bond strength is kq ± Ak with equal probability. 
This is motivated from a theorem 0] which states that 
in two dimensions any amount of bond disorder will be 
sufficient to soften a first order transition into a second 
order transition. Secondly, we have studied a model were 
only half of the sites have next-nearest neighbour inter- 
action. This model will clearly share many of the quali- 
tative features of the original model, but the effect of the 
next-nearest neighbour interactions will be reduced. 



II. QMC SIMULATIONS 

We have performed Quantum Monte Carlo simula- 
tions using the Stochastic Series Expansion fSSE)p^Lll^| 
method. In the SSE method the Hamilton operator is 
written as a sum of bond operators 



(3) 



The sum in Eq.[3|is over all the bonds on the lattice, Hd,b 
is an operator working on bond 6, which is diagonal in 
the basis chosen to represent the spin space, and H Q d t b 
is an off-diagonal operator. For spin models with z axis 
magnetization as basis, the operator Hd,b will be 



Hd,b — J SnySjQ,), 



(4) 



where i(b) and j(b) are the two sites connected by bond 
b. H a d,b is an off-diagonal operator, and in the case of 
spin models we will have H Q d t b given by 



H od ,b - - ( s t(b) S 3 (b) + S i(b) S j~(b) 



(5) 



Observe that for the actual simulations the operators 
Hd,b are scaled and shifted |T3| to ensure 



The formal expression for the partition function is then 
expanded, which yields the following representation 



{a} n {S„} 



IK 



(7) 



Here S n is a sequence of n pairs, each pair consisting of 
a variable denoting operator type and a bond index, i.e. 



S n = { (ai,h), (a 2 ,b 2 ) , ... (a n ,b n ) }. 



(8) 



The variable a, in Eq.[Sl denotes type of operator and can 
be either diagonal or off-diagonal. The SSE method then 
consists of doing importance sampling in the combined 
space |{a}) <g>5 n . The actual updates are of two different 
types. The diagonal updates insert or remove a diago- 
nal operator -ff^b, thereby changing the expansion order 
n — ► n ± 1. The off-diagonal operators change operator 
types Hd.b <-> H dfi and flip the corresponding spins, this 
must be done in a way which ensures periodicity in the 
(3 direction, i.e. |a(0)) — \a(n)). For the off diagonal 
updates, the advent of loop updates 01 nas si gnifi cantly 
improved the performance of SSE simulations |l8lll^| . 

For the ordinary 5 = 1/2 Heisenberg model, SSE simu- 
lations with operator loop update are particularly simple. 
In order to include the next-nearest neighbour interac- 
tions, we need to modify the algorithm slightly. In the 
case of the diagonal updates this merely amounts to in- 
cluding the extra factor k in the weight calculation for 
the next-nearest neighbour bonds. Whereas for the op- 
erator loop the next-nearest neighbour interactions have 
a more profound effect. These interactions are only diag- 
onal, i.e. the incoming and outgoing spin states must be 
equal. Furthermore, the next-nearest neighbour bonds 
can only connect antiparallel (k > 0) spins. The re- 
sult of this is that the next-nearest neighbour bonds 
"freeze" a substantial part of the spin configuration, and 
only those spins/operators not directly linked to a next- 
nearest neighbour bond are amenable to operator loop 
update, as illustrated in Fig. ^ Clearly, this freezing 
affects the performance of the simulations in a negative 
way, in particular for intermediate values of k. 



A. Observables 

To differentiate between the different types of order in 
the model, we have studied the structure factor 



S(Q) = ^((j2(a\S z (r)\a)e^ 



HdA li) = | U) 



H d ,b\ TT) = 0. (6) 



(9) 



for different values of Q. An estimator of 5(Q) taking 
all the intermediate SSE states into account can be found 
in [16j . For the remaining part of the text we will make 
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FIG. 1: A chain of six spins, depicted with an operator se- 
quence of length n = 8 in the direction. The filled bars 
denote diagonal interactions SfSj and the open bars are flip 
operators S^~S~ + S~ . The dashed region shows spins/p- 
slices which have been frozen by the next-nearest neighbour 
bond in the middle. 

frequent use of the terms staggered and striped magne- 
tization, these quantities are defined as 



= (10) 

Mf W)0) - jjy/S(Tr,0) + S(0,n). (11) 

The upper index indicates that the magnetization is eval- 
uated along the z axis, and the lower index is the direc- 
tion of Q in the evalution of Eq.EI i.e. (it, ir) for staggered 
and (tt, 0) for striped magnetization. The factor of three 
in Eg. 1101 is included to account for rotational averaging 
among the three directions in spin space. When k is fi- 
nite, isotropy in spin space is explicitly broken. We have 
nevertheless retained this factor to get continuous formu- 
lae around k = 0. In addition to the structure factor, we 
have also measured two other quantities, namely specific 
heat and superfluid density. 
The specific heat is given by 

Cv=P 2 ^\nZ={n(n-l))-(n) 2 . (12) 

Here, n refers to the summation variable in Eq. This 
summation is truncated in a stochastic manner and n 
is thus promoted to a dynamical variable in SSE. Wc 



will not exhibit results for Cy explicitly, but have used 
the anomalies in this quantity to corroborate the phase 
boundaries shown in Fig. 0] (with the exception of the 
line separating the superfluid phase from the disordered 
phase, see comments on this below). 

The estimator for the superfluid density (0(2) order- 
ing) is given by [16| 

ps = m - ^ )2) + < w - ^ )2) ) ' (13) 

where N+/N~ is the number of S~ and S~S^ oper- 
ators applied along bonds in the \i direction. 

III. RESULTS 

As mentioned previously, the model in the form of 
Eq. |21 has already been studied in Ref. . As a bench- 
mark of our QMC methods, we started with this model 
to reproduce the results of Ref. Figs. [3 and |3 shows 
the staggered and striped magnetization along the z axis, 
as a function of k. 
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FIG. 2: (Color online) The staggered magnetization along the 
z axis in the ground state ((3 — 10), as a function of k. 

From these two figures we conclude the following. (1) 
For negative k the magnetization along the z axis is en- 
hanced; this is expected since the k < system is not 
frustrated. (2) For k > the magnetization is immedi- 
ately tilted away from the z axis, leaving zero magnetiza- 
tion along the z axis. (3) For large k the magnetization 
again orders along the z axis, in this case in a striped 
formation. The transition to the striped phase is a first 
order transition, the discontinuous jumps in M z (tt, 0) in 
Fig. |3 indicate this, and it is also confirmed by a more 
detailed analysis of histograms of e.g. the striped or- 
der parameter or number of next-nearest neighbour 
bond operators. Hence, for this case the order-to-order 
transition depicted in Fig. (a transition from super- 
fluid order, equivalently 0(2) order, to stripe order) falls 
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FIG. 3: (Color online) The striped magnetization in the 
ground state (/3 = 10), as a function of n. There is a first 
order transition at re ~ 1.205. 



within the standard Landau- Ginzburg- Wilson paradigm 
of phase transitions. 

For k < and k > the ordered state breaks a discrete 
symmetry, and the order persists for finite T. In the in- 
termediate regime, < k < 1.205 the remaining model 
is an antiferromagnetic 2DXY model with a continuous 
0(2) symmetry. According to the Mermin Wagner theo- 
rem this symmetry can not be spontaneously broken at 
finite T. However, there is finite spin stiffness and topo- 
logical order up to a temperature Tbkt where the order 
vanishes in a Berezinski-Kosterlitz-Thouless transition. 
The boundary of this region has been (approximately) 
located by equating ps{T) with 2T/tt. All in all we have 
found the phase diagram presented in Fig. 0] for Eq. 



0.6 



0.4 



0.2 





Disordered phase 




Neel order \ 




Striped order 




Superfluid 





-0.5 



0.5 



1.5 



FIG. 4: Phase diagram for the model in the k, T plane. For 
k > k c , and k < the model has magnetic order, striped and 
staggered respectively. In the intermediate /t range values 
there is no finite T magnetic order, however there is a super- 
fluid order which persists into the finite T region. The line 
separating topological 2DXY order from the normal phase 
is determined with considerably less precision than the two 
other lines. 



A. Disordered system 

The effect of disorder on phase transitions is a much 
studied topic. In the case of continuous transitions, the 
Harris criterion pol l2l| states that disorder will change 
the universality class of the transition, i.e. be relevant, if 
the exponents of the pure system satisfy v < 2/d. In 
the case of first order transitions, disorder can soften 
the transition into a continuous transition, in the case 
of two dimensions any amount of disorder is sufficient 
|14|. whereas a finite amount is needed in three dimen- 
sions. These predictions have been confirmed for the 
Potts model in both two and three dimension]^, l23| . 
We have, however, not found tests of these predictions 
for a first order quantum phase transition. We have in- 
vestigated what happens with the first order quantum 
phase transition at K rs 1.205 when disorder is included 
in the model. Along each bond is k is given by 



K = Kq ± Ak, 



(14) 



with equal probability. We have focused on the striped 
order parameter M z (tt, 0) in the vicinity of k ~ 1.205, in 
order to compare with Fig. [21 In the disordered system 
we must perform both ordinary thermodynamic averag- 
ing and subsequently disorder averaging. In e.g. Fig. [S] 
the plotted quantity is given by 



1 



N 



M*(7T,0) = -^Mf(7r,0), 



(15) 



where M?(tt, 0) is the striped magnetization in disorder 
realization i, calculated according to Eq. ^] and N is 
the total number of disorder realizations. The number of 
disorder realizations has typically been N = 100. 

Fig. El shows the disorder averaged striped magneti- 
zation as a function of no for Ak = 0.05. The strongly 
discontinuous features of M z (tt, 0) from Fig.|3]are washed 
out when disorder is introduced. From this, we conclude 
that the transition changes order when disorder is in- 
troduced. The location of the critical point coincides 
with the original transition point of the uniform system. 
We have not varied Ak systematically, our results (not 
shown) indicate that the system is not very sensitive to 
Ak variations. 

The low-K region in Fig. [S] is an 0(2) ordered state. 
The large-K region is an Ising-ordered state with an addi- 
tional stripe order. Hence, this quantum phase transition 
is a transition from an ordered state to another ordered 
state, and since it is continuous, it represents a candidate 
example of so-called deconfined quantum criticality. 



B. AB system 

For large k the Ising next-nearest neighbour interac- 
tion in Eq. [21 is a very strong interaction. In an attempt 
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FIG. 5: (Color online) Disorder averaged value of M z (n, 0) as 
a function of Ko for two different system sizes. The indicated 
error bars are standard error estimates from the independent 
disorder realizations. The temperature coupling is ft = 10. 



where the sum is only over A sites. Because the sum is 
limited to the A sites, full polarization corresponds to 
M z = 0.25. Fig. shows the staggered magnetization 
among the A sites as a function of k. 
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to soften the transition to the high k state into a con- 
tinuous transition we have devised a model consisting of 
two "atom" types A and B, where the Ising next-nearest 
neighbour interaction is only between the A atoms, the 
scenario is illustrated in Fig. EJ 




FIG. 6: The binary system consisting of two different atomic 
species. Only type A (•) has next-nearest neighbour interac- 
tion, illustrated with diagonal lines. 

This model has many of the same qualitative proper- 
ties as the original model, in particular small values of k 
will frustrate the system and tilt the magnetization down 
in the xy plane. For large k the A atoms will form an 
antiferromagnetically ordered state, with AF magnetiza- 
tion along the z axis. In this state the A and B sites 
will decouple, and the B sites will be disordered with no 
net contribution to the energy of the system. As an or- 
der parameter for this transition we have considered the 
staggered magnetization along the z axis, for the A sites, 
i.e. 



M&wn = WSA(7r/2,7r/2), (16) 



FIG. 7: (Color online) The staggered magnetization of the A 
sites in the groundstate of the AB model, as a function of k. 
The solid line is a least squares fit to a power law. 

Comparing with the striped magnetization of the uni- 
form model, Fig. |3 we see that the critical coupling 
k c w 1.705 of the AB system is much larger than for 
the uniform system. This is reasonable, since the next- 
nearest neighbour interaction only operates on half of the 
sites. Furthermore, Fig. [7| shows that the transition to 
the ordered state is much smoother than in the uniform 
model, hence the figure indicates that the transition is 
continuous. The continuous nature of the transition is 
also confirmed by considering the number distribution of 
e.g. next-nearest neighbour operators at critical point. 
This quantity is unimodal, whereas for the uniform model 
it is strongly bimodal, reminiscent of a first order tran- 
sition. From this, we conclude that the AB-model defor- 
mation of the Hamiltonian in Eq. [21 suffices to promote 
quantum deconfined criticality. Stripe order is lost by re- 
moving every second next-nearest neighbor coupling, but 
an order-to-order quantum phase transition nevertheless 
remains. Namely, the transition is that from an 0(2) (su- 
pcrfluid) ordered state at low and intermediate values of 
ft to a Zi (Ising) ordered state at high values of n. Since 
the transition is second order, it represents a second can- 
didate example of deconfined quantum criticality. 

We have not made attempts at completely determining 
the critical exponents at the transition. However a fit of 
M z from the L = 32 system to the functional form 

M (W/2) « l«- K c| /3 (17) 

with k c = 1.705 gave good results, with a critical expo- 
nent (3 « 0.25. The fit is shown as a solid line in Fig. [7J 
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IV. CONCLUSION 

We have studied varieties of the antiferromagnetic 
Heisenberg model in two dimensions with additional 
next-nearest neighbour Ising exchange. This model has a 
strong first order transition at n c m 1.205. We have stud- 
ied two generalisations of the model, one based on dis- 
ordered next-nearest neighbour couplings, and another 
where only half the sites are endowed with next-nearest 
neighbour interaction. Both the generalised models fea- 
ture continuous quantum phase transitions from one or- 
dered state to another ordered state. As such, these two 
cases represent candidate examples of deconfined quan- 
tum criticality. Frustrated interactions is an essential 



part of the models we have considered, and as such they 
are distinct from the model already considered by Sand- 
vik et.al. 00; where (possible) deconfined criticality is 
brought about by ring-exchange. 
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